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A theory of quantum-mechanical generation of cosmological perturba- 
; tions is considered. The conclusion of this study is that if the large-angular- 

scale anisotropy in the cosmic microwave background radiation is caused 
by the long-wavelength cosmological perturbations of quantum mechanical 



origin, they are, most likely, gravitational waves, rather than density per- 
turbations or rotational perturbations. Some disagreements with previous 
publications are clarified. 



This contribution to the Proceedings is based on Reference [34]. 

Sh '. 

(3JT), 1. Introduction 

The ongoing and planned high-precision measurements of the anisotropy in 
^ the cosmic microwave background radiation (CMBR) [1] may have serious 

impact on our views about the very early Universe. The large-angular-scale 
anisotropy in the CMBR is most likely caused by cosmological perturba- 
tions with wavelengths of the order and longer than the present day Hubble 
radius Ih- It is reasonable to expect that so long- wavelength perturbations 
are "primordial", survived from the epochs when the Universe was much 
younger. The wavelengths of the perturbations have enormously grown up 
since the time of generation but other physical characteristics of the pertur- 
bations can still carry imprints of their origin. This cannot be said with the 
same degree of certainty about the relatively short-wavelength cosmolog- 
ical perturbations (unless they are gravitational waves) which could have 
been distorted and contaminated in course of their life by many physical 
processes occurring in the Universe. 
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It is remarkable that the origin of all three possible types of cosmological 
perturbations, that is the origin of density perturbations, rotational per- 
turbations, and gravitational waves, may be of purely quantum-mechanical 
nature. Cosmological perturbations can be treated as excitations in grav- 
itational field. In case of gravitational waves, they are just excitations in 
gravitational field itself. In case of density perturbations and rotational 
perturbations, they are excitations in gravitational field which accompany 
excitations in matter. In the very distant past, the density and rotational 
perturbations were excitations in the primeval medium that was filling the 
Universe at that time. The quantum-mechanical generation mechanism of 
cosmological perturbations relies only upon the existence of their zero-point 
quantum fluctuations and the nonvanishing parametric coupling of the per- 
turbations to the variable gravitational field of the homogeneous isotropic 
Universe. The strong variable gravitational field of the very early Universe 
played the role of the pump field. It supplied energy to the zero-point quan- 
tum fluctuations and amplified them. More precisely, the initial vacuum 
quantum state of each mode of the perturbations has been transformed, 
as a result of the quantum mechanical Schrodinger evolution, into a multi- 
particle state known as a squeezed vacuum quantum state. The generated 
perturbations have formed a collection of standing waves. The gravitational 
field of each of the three types of quantum mechanically generated per- 
turbations can affect the propagating photons of the CMBR and produce 
anisotropy in CMBR. 

It was already emphasized [2] that there is a significant qualitative dif- 
ference between gravitational waves on one side and density and rotational 
perturbations on the other side, with regard to possibility of their quan- 
tum mechanical generation. Gravitational waves oscillate in the absence of 
external gravitational fields, and their appropriate parametric coupling to 
the pump field follows directly from the Einstein equations. The parametric 
excitation vanishes only if the cosmological scale factor obeys the equation 
a" /a = 0, that is when there is no any pump field at all, a(rj) = const, or 
when the coupling a' is time independent. Up to this exception, one can say 
that the quantum mechanical generation of gravitational waves (relic gravi- 
tons) is unavoidable [3]. As for density and rotational perturbations, they 
are perturbations in matter being accompanied by perturbations of gravi- 
tational field. The ability to support oscillations of density and/or rotation 
and the form of their coupling to the pump field depend on a particular 
model of matter and its energy-momentum tensor. The very possibility of 
the quantum mechanical generation of these perturbations is model de- 
pendent. Recalling Einstein's definition of two pillars supporting general 
relativity, one can say that the quantum mechanically generated gravita- 
tional waves are associated with the pillar made of marble, while density 
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and rotational perturbations are associated with the other one. 

Presently, we do not know if the observed large-scale structure in the 
Universe did really develop from the primordial density perturbations of 
quantum-mechanical origin. It is conceivable that the primordial pertur- 
bations arose as a result of, say, a sophisticated nonsimultaneous phase 
transition in the very early Universe or as a consequence of the hypotheti- 
cal cosmic defects. All reasonable hypotheses should be carefully studied. It 
is necessary, however, to see clearly the difference in the number and value 
of initial assumptions: the quantum- mechanical generating mechanism re- 
quires only the validity of general relativity and basic principles of quantum 
field theory. Moreover, the analogous process of parametric amplification 
and generation of squeezed vacuum states of light has been demonstrated 
in quantum optics laboratory experiments. It seems that we are given a 
great opportunity to derive the cosmological perturbations from the very 
first principles. This is why the quantum-mechanical generation mechanism 
will be the subject of our study here. 

As was discussed above, gravitational waves are generated universally in 
contrast to density and rotational perturbations which require appropriate 
matter. A particular sort of matter that has received much attention in the 
recent cosmological literature, especially in the literature on inflation [4], 
is one or another version of a scalar field. Scalar fields is a nice theoretical 
model that has been used in physics in many different studies. Whether 
the global scalar fields do really exist in nature and, if so, whether they 
couple to gravity in the way we want, is presently unknown. However, we 
will follow the modern tradition in theoretical physics which states that 
everything that is not forbidden is allowed. One can at least guarantee that 
a sort of inflationary expansion is a typical feature (attracting separatrix) 
in the space of homogeneous isotropic solutions to the Einstein equations 
with certain scalar fields [5] . Scalar fields cannot support rotational pertur- 
bations but they can support density perturbations. 

Specifically, we will study a scalar field ip(t, x 1 , x 2 , x 3 ) with the energy- 
momentum tensor 



(1) 



where V(tp) is an arbitrary scalar field potential, comma denotes a par- 
tial derivative. This model of the primeval cosmological medium satisfies 
both conditions for the quantum mechanical generation of density pertur- 
bations be possible. First, the field can obviously support free oscillations in 
Minkowski space-time. Second, the explicit form of the energy-momentum 
tensor (1) reflects the appropriate (minimal, the same as for gravitational 
waves) coupling of the scalar field to gravity which was chosen by our will. 
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So, on general grounds and by analogy with gravitational waves, one can ex- 
pect that some amount of density perturbations might have been generated 
by strong variable gravitational field of the early Universe. The problem is 
to quantify this expectation and to derive the observational predictions, as 
reliable and detailed as possible, including the expected variations in the 
CMBR. 

Scalar fields and scalar field perturbations is a very popular subject in 
the framework of inflationary cosmologies. So popular, that many believe 
that the inflationary type of expansion is conditioned by the existence of 
scalar fields and that the very possibility to generate perturbations quan- 
tum mechanically relies on the existence of the De Sitter event horizon. 
This is not so. Inflation, if understood as a statement about the behavior 
of the time dependent cosmological scale factor, and not about creating 
and resolving the particle physics paradoxes, is a phenomenon more gen- 
eral than one particular realization of it with the help of a scalar field. [The 
attitude toward the archetype inflationary solution — exponential expan- 
sion — has changed over the years. Astronomers of the older generation 
were embarrassed with the De Sitter solution but tried to apply it for the 
explanation of the galaxies' red shifts and statistics of quasars in the most 
recent Universe. Cosmologists of our time take the exponential expansion 
as something almost proven but apply it to the very remote stages of evo- 
lution, somewhere near the Planck time.] And the quantum mechanical 
generation of perturbations is a phenomenon more general and universal 
than such concepts as global scalar fields, event horizons, and inflation. 
If it turns out that the inflationary hypothesis contradicts observations, 
the quantum mechanical generating mechanism will not die together with 
inflation. There is little doubt, for instance, that the search for relic gravita- 
tional waves will continue, with may be larger emphasis on relatively short 
waves rather than on long waves [6] . And a test of the quantum mechanical 
origin of cosmological perturbations will be a test of their origin, not a test 
of inflation specifically. 

The generation of density perturbations in inflationary models governed 
by the scalar field (1) was a subject of discussion in many research and re- 
view articles, and books. If one consults the most recent literature, one 
can find that the current situation is often summarized in the following, or 
similar, words: "Exponential inflation predicts a scale-invariant, Gaussian 
spectrum of scalar fluctuations and a smaller amount of tensor fluctu- 
ations ... . Other inflationary models, for instance power-law inflation 
predict spectra slightly tilted away from scale invariance." (See, for ex- 
ample, [7] and references therein.) The expected amplitudes of density 
perturbations are usually quoted in the following, or equivalent, form (see, 
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for instance, [8] and references therein): 

(Sp\ hm _ mn 2 H 2 (<p) 
Wa ~ 87r3/2 \ H '&)\ ' 

where the quantities on the right-hand-side are supposed to be evaluated 
"when the scale A crossed the Hubble radius during inflation" . The denom- 
inator of this expression depends on the derivative of the Hubble parameter 
and goes to zero in the limit of exponential inflation. Apparently, this for- 
mula says that the predicted amplitudes of the scale-invariant spectrum 
are arbitrarily close to infinity, and the amplitudes of nearby spectra are 
"slightly tilted away" from infinity. According to this formula, the ampli- 
tudes of density perturbations are many orders of magnitude larger than the 
amplitudes of gravitational waves, if the expansion rate is sufficiently close 
to the exponential inflation. The belief that the amplitudes of density per- 
turbations are larger, or much larger, than the amplitudes of gravitational 
waves is considered to be a strong prediction of inflationary models based 
on the scalar field (1). For instance, the author of Ref. [8] concludes: "An 
observation violating this condition at any scale would immediately rule 
out the general class of models we are considering". The expected contri- 
bution of density perturbations and gravitational waves to the quadrupole 
anisotropy of CMBR was also under study. The authors of Ref. [9] (see also 
[10] and references therein) say that "The ratio of gravitational wave (T) 
to energy-density perturbations (S) contributions to the CMB quadrupole 
anisotropy is predicted to be T/S = 21(1 +7)", where 7, in that paper, is 
the parameter in the equation of state for matter governing the inflationary 
expansion, p = jp. According to this formula, T/S vanishes in the limit 
of 7 = —1, that is in the limit of strictly exponential (De Sitter) inflation. 
Apparently, this formula for T/S is based on the authors' assumption that 
the effectiveness of generation of density perturbations is the higher the 
closer the expansion law to the exponential inflation, and goes to infinity 
in the limit of 7 = —1. Apparently, this is why T/S goes to zero in this 
limit. The statements about density perturbations are sometimes charac- 
terized as such that have been "widely studied and there is broad agreement 
regarding both methods and results ..." [11]. 

I suspect that the present work will not belong to that category of stud- 
ies that enjoyed the "broad agreement"; my conclusions are considerably 
different from what was described in the preceeding paragraph. I will be 
arguing that there is no linear density perturbations at all at the purely 
exponential (De Sitter) inflationary stage for models governed by the scalar 
field (1). Density perturbations can only arise as a result of violation of the 
purely exponential expansion and transition to the radiation-dominated 
stage. Regardless of how close to zero was the derivative of the Hubble 
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parameter "when the scale A crossed the Hubble radius during inflation" , 
the today's amplitudes of density perturbations are finite. The amplitudes 
of gravitational waves are typically a little larger than the amplitudes of 
density perturbations, at least in the long wavelength limit where spectra 
are smooth and have the power-law behavior. Correspondingly, the contri- 
bution of density perturbations to the quadrupole anisotropy is never much 
larger than the contribution of gravitational waves. In fact, it is somewhat 
smaller in the limit of long waves. 

As for the statistical properties of cosmological perturbations and, hence, 
the statistical properties of the CMBR fluctuations caused by them, it was 
already emphasized [12] that they are determined by the statistics of quan- 
tum states being generated, namely by the statistics of squeezed vacuum 
quantum states. 

In this paper, we will present detailed derivations of equations and their 
solutions in order to make it possible for the interested reader to compare 
the present calculations with those of other authors. The structure of the 
paper is the following. In Sec. II we present the general equations for cos- 
mological perturbations. The equations are applicable for matter with the 
energy-momentum tensor of arbitrary form. They use only the defining 
property of cosmological perturbations, namely that the perturbations arc 
based on the scalar, vector, and tensor functions of spatial coordinates. In 
Sec. Ill we apply these equations specifically to the initial stage (i stage) 
of cosmological evolution which is assumed to be governed by the scalar 
field with the energy-momentum tensor (1). The main purpose is to study 
the density perturbations. No assumptions about a particular form of the 
scalar field potential V(<p) or a particular (for example, inflationary) type 
of expansion are being made a priori. The time-dependent coefficients of 
the differential equations for the perturbations are expressed in terms of the 
scale factor (and its derivatives) only. This reflects the underlying interac- 
tion of the perturbations with the variable gravitational pump field. The 
determination of all unknown functions describing density perturbations is 
reduced to solving a single differential equation which is very similar to the 
equation for gravitational waves. The behavior of solutions during a more 
or less gradual transition from the i stage to the radiation-dominated stage 
(e stage) is studied in Sec. IV. In order to deal with simple exact solutions 
at both stages we will be interested in a sharp transition from the i stage 
to the e stage. In Sec. V we apply the perturbation equations to the perfect 
fluid matter with arbitrary velocity of sound. We present solutions to these 
equations at the e stage and the matter-dominated stage (m stage) in the 
form convenient for matching the solutions at all three stages (i, e, m). In 
Sec. VI we join the solutions, find the coefficients which were undetermined 
so far, and express the solution at the m stage entirely in terms of the func- 
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tions (and their first time derivatives) describing the perturbations at the 
time of joining the i and e stages. As a preparation for quantization of the 
perturbations, we briefly discuss the density and rotational perturbations of 
matter placed in the Minkowski space-time, that is neglecting gravitational 
fields (Sec. VII). The quantization of density perturbations is performed 
in Sec. VIII. This procedure essentially repeats the steps which have been 
previously done for gravitational waves and rotational perturbations [12,2]. 
The quantum mechanically generated perturbations are placed in squeezed 
vacuum quantum states. Classically, one can think of the perturbations as 
of a stochastic collection of standing waves. The justification and necessity 
of the so-called Sakharov's oscillations in the spectra of density perturba- 
tions is discussed. In order to get the analytic results as detailed as possible, 
we specialize the scale factor of the i stage to the ry-time power laws which 
include inflationary models. This allows us to derive concrete power-law 
spectra of density perturbations at the m-stage. In Sec. IX we derive an ex- 
act formula for the angular correlation function of the CMBR temperature 
variations 5T/T caused by squeezed density perturbations. The multipole 
decomposition of the correlation function begins from the monopole term. 
The contributions to the monopole and dipole terms produced by individ- 
ual waves with wavelengths exceeding Ih are strongly suppressed, which is 
in agreement with previous results [13]. Nevertheless, one should be aware 
that not only the entire quadrupole but also some little portions of the 
measured mean temperature of CMBR and its dipole variation may be 
caused by perturbations of quantum mechanical origin. For one and the 
same cosmological model, the contribution of density perturbations to the 
quadrupole anisotropy of CMBR is never much larger than the contribution 
of gravitational waves, at least for models considered here. 

2. General Equations for Cosmological Perturbations 

The unperturbed spatially-flat FLRW (Friedmann-Lemaitre-Robertson- 
Walker) cosmological models are described by the metric 

ds 2 = -c 2 dt 2 +a 2 (t)(dx l2 +dx 22 +dx 32 ) = -a 2 (r))(dr) 2 -dx l2 -dx 22 -dx 32 ) . 

(2) 

Here we restrict ourselves to study of the spatially-flat models. This is 
done for purposes of simplicity only, not because of any inflationary prej- 
udices. If it turns out that the density parameter Q is, say, Q = 0.8, in- 
stead of Q = 1, the inflationary hypothesis will prove to be wrong but the 
quantum-mechanical generating mechanism will still be perfectly well effec- 
tive, although with certain modifications as compared to the spatially-flat 
case. 
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The scale factor a(rj) is governed by matter with the unperturbed values 
of energy density eo, T° = — eo, and pressure po, if = poSf: 




KC 



= Wo (3) 

U~ \ U I \ U J 

where k = 87rG/c 4 and a prime is d/drj, d/drj = (a/c)d/dt. The Hubble 
parameter is H = a /a = ca' /a 2 . 

It is convenient to introduce two new functions of the scale factor: 

a' a' 
<x(ri) = - , 7(»7) = 1 - — • ( 4 ) 
a a z 

In terms of t-time the function 7 is j(t) = —(H/H 2 ). Due to Eqs. (3) one 
has 

K(e +po) = —5-7 • (5) 
a 1 

The function 7(77) becomes a constant if a{rf) is governed by matter with 
the effective equation of state po = qeo, where q = const. The scale factor 
a(rj) takes on the ry-time power law behavior 

a = l \v\^ (6) 

(77 must be negative for expanding models with l+(5 < 0) where l a and (3 are 
constants. (This parameterization has been used before [3].) The constant 
l Q has the dimensionality of length. It follows from Eqs. (3) and (4) that 
7 = (2+/3)/(l+/3), q((3) = (1— (3)/3(l+(3), where the parameter (3 can vary in 
the interval — 00 < [3 < 00. In particular, 7 = 2 at the radiation-dominated 
stage, and 7 = 3/2 at the matter-dominated stage. Note that 7 = in 
case of purely exponential (De Sitter) expansion for which a(t) ~ e Ht , 
H = const, a(r]) = / |r/| -1 , f3 = —2. 

One can also derive from Eqs. (3) the relationship 

f- 1+ r-sf->^- (7) 

The ratio p'o/e' becomes a constant for the scale factors (6), namely: 
Pq/cq = q{0). In particular, p' /e' goes to —1 in the limit of (3 = —2. 

Let us now turn to the perturbations of the spatially-flat FLRW models. 
Without resticting in any way the physical content of the problem, it is 
convenient to work in the class of synchronous coordinate systems. (At this 



9 



point the reader may have to be ready to exhibit certain resistance to the 
pressure from the proponents of the "gauge-invariant" formalisms.) Using 
the r/-time coordinate, one can write the general expression for the metric 
tensor including perturbations as: 



ds = —a 



drj 2 — (5{j + hij)dx L dx^ 



where the six functions hij describe all the independent components of the 
perturbed gravitational field. 

The construction of cosmological perturbations [14,15] (see also [16]) 
is based on the classification of the eigenfunctions to the 3-dimensional 
Laplace operator. Density perturbations are associated with the scalar func- 
tions Qix 1 , x 2 , x 3 ) satisfying the equation 

Q>\ t+ n 2 Q = (8) 

valid in three-space dl 2 = dx l2 + dx 2,2 + dx 3,2 . For each wave vector n, one 
can choose two linearly independent solutions to Eq. (8) in the form e mx 
and e~ mx . From a given scalar field Q one can construct a vector field 
and two tensor fields: SikQ and Q t i t k = —nin^Q. For each n, the general 
perturbation of the energy-momentum tensor and the accompanying per- 
turbation of the gravitational field can be written as a sum of products of 
time-dependent amplitudes and spatial functions introduced above. 

The components of the perturbed gravitational field can be written as 

= hQ5ij + hin~ 2 Q ti j (9) 



where the function h(rj) represents the scalar (proportional to <%Q) per- 
turbation of the gravitational field while the function hi (rj) represents the 
longitudinal-longitudinal (proportional to niUjQ) perturbation. The gen- 
eral expression for including perturbations can be written as 

1 1 , 1 

T° = -€q n ei( 2> T i=^2^,Q,i^ T o = o^Q' 1 ' 

a z a z a z 

T t h = PoSf + ^{pi+Pi)Qd^ + ^n- 2 Pl Q ti ' k . (10) 



The form of Eqs. (9) and (10) is based solely on the definition of density 
perturbations and our choice of synchronous coordinate systems. In all 
other respects, the representation (9) and (10) is general. The particular 
notations for arbitrary functions describing the perturbations are chosen 
for later convenience. 
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The arbitrary functions h(rj), hi(rj), ei(rj), p\{jfj Pi(r)), C'(v) should sat- 
isfy all together the perturbed Einstein equations: 

3ah' + n 2 h — ah\ = Ke± (11) 
ti =kC (12) 
—h" — 2ah' = npi (13) 

-(h'{ + 2ah' l -n 2 h) = k Pi . (14) 

There are too many unknown functions to be found from Eqs. (11)-(14). 
This requires us to specify a model for matter and its energy-momentum 
tensor. We will consider three consecutive stages of expansion: i stage 
governed by the scalar field (1), and the subsequent e and m stages governed 
by perfect fluid with the energy-momentum tensor 

V = + ( 15 ) 

The equation of state at the e and m stages is p = |e and p = 0, respectively. 

Now we will turn to the rotational perturbations and gravitational 
waves. The rotational perturbations are based on the vector functions Qi 
satisfying the equations 

Qi,k' k + n 2 Qi = Q\ t = ■ 

The nonvanishing components of the perturbed gravitational field and the 
nonvanishing components of the perturbed energy-momentum tensor can 
be written in the form 

hij = 1h(ji)A l3 , T° = -Tl = a~ 2 co( V )Qi , T t k = 2 a - 2 n - l X { 11 )A k % 

where A; L j = — ^(Qij + Qj,i)- The arbitrary functions h(rf), u(rf), x(v) 
should satisfy together the perturbed Einstein equations: 

h" + 2ati = 2ku- 1 x{v) , 
~h' = kuj(ti) . 

And, finally, gravitational waves are based on the tensor functions Qij 
satisfying the equations 

Qij,k' k + n 2 Qij = , Q ij j = 0, Q ll 5 t3 = . 

The components of the perturbed gravitational field should obey the per- 
turbed Einstein equations 



hi" + 2ah{ + n 2 hi j = 
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Note, that the components Tj fe of the perturbed energy-momentum tensor 
vanish for simple models of matter, such as perfect fluids, but, in general, 
they may be different from zero. However, we put them equal to zero since 
we are not interested in the generation of gravitational waves by anisotropic 
stresses of the material medium. Instead, we are interested in the quantum- 
mechanical generation of gravitational waves by variable pump field. This 
process is independent of the material anisotropic stresses and takes place 
even if there is no classical sources at all. 
By introducing 

hij = - ^{v)Qij 

one can reduce the basic equation to the form of the equation for the 
parametrically excited oscillator with the time-dependent frequency 

which reflects the underlying physics of the phenomenon. This equation has 
been studied before (for a review, see Ref. [19]). 

At the first glance, the original equations for rotational and density 
perturbations do not look like equations for a parametrically coupled oscil- 
lator. However, they can be reduced to this equation for appropriate models 
of matter whenever the amplification process is at all possible. In case of 
rotational perturbations this has been done in Ref. [2]. For density pertur- 
bations we will do it below by concentrating on the scalar field matter with 
the energy- momentum tensor (1). 



3. Density Perturbations at the Initial Stage of Expansion 
Governed by a Scalar Field 

At the i stage, the evolution of the scale factor a(rj) is determined by the 
unperturbed homogeneous scalar field ip = ip (v)- The unperturbed values 
eo, Po are given by Eq. (1): 

60 = i^o) 2 + ^M (16) 

Po = ^(Vo)-^) • (17) 

By summing up Eqs. (16) and (17) and comparing the result with Eq. (5) 
one can derive the equation 



K (<p' ) 2 = 2a 2 7 



(18) 
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It follows from this equation that 7 > for the scale factors governed by 
the scalar field (1). The De Sitter case corresponds to <p' = 0, ipo = const 
and eo = — Po = V(<po) = const. 

If <p' 7^ 0, one can use the equation 

e' = -3a(e + Po) 

which is a consequence of Eq. (3), and obtain with the help of Eqs. (16) 
and (17): 

¥>{,' + 2av' + a 2 V v = (19) 

where = dV(ip)/dip, the derivative is taken at 99 = ifQ. The further 
useful relationships following from Eqs. (18) and (19) are: 

4 = U\± (20) 

-4^ = ^ + ^ + . (21) 

(f/ a 2 7 

The perturbations of the gravitational field are associated with the per- 
turbations of the scalar field. We will write the perturbations of the scalar 
field as 

<p = Mv) + Mv)Q • (22) 

Having at our disposition the energy-momentum tensor (1) and the def- 
initions (22), (9), (10) we can directly calculate the functions ei, pi, 
Pi- 

ei = (p' ip[ + a 2 (fiV^ (23) 
£' = (24) 

Pi = <Po<Pi - «Vi^ (25) 
Pi = . (26) 

We will now assume that <^q 7= 0. The De Sitter case (p' = will be 
considered separately at the end of this Section. It follows from Eq. (24) 
that ip\ = — £'/Vo- Inserting this value of 991 into Eqs. (23) and (25) and 
using Eqs. (20) and (21), one can express e\, pi, in terms of h{rj): 



Kei = -h" + ti\2a + 2^- + j-j (27) 
Kpi = -h"-2ati . (28) 
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We should now return to the perturbed Einstein equations (11)-(14) 
making use of Eqs. (27) and (28). Equation (11) can be written as an 
expression for h[{ri) in terms of h(rj): 



H = T 



a 



h" l/i a -2--^)+n 2 h 

a 7 



(29) 



Equation (12) expresses £' in terms of h'. Equation (13) is satisfied identi- 
cally. Equation (14) reads as 



h" + 2ah[ -n 2 h = 



(30) 



Thus, if one knows the function h{rj), all other functions describing the 
density perturbations, namely hi(r]), Pi(rj), £1(77), and £'(77), can be found 
with the help of Eqs. (29), (28), (27) (or, equivalently, (11)) and (12). To 
derive the equation for 71(77) we substitute Eq. (29) into Eq. (30) and obtain 



ti"+h"[3a>y-- \+ti 



7 



n 2 -2a'+2 7 a 2 -^- - 
07 V 7 



+n 2 a-fh = 0. (31) 



Equation (31) is a third-order differential equation. There should be no 
wonder (and no panic) on this occasion. One of solutions to this equation 
we know in advance, this solution is 



h = C™ 
a 



(32) 



where C is an arbitrary constant. We could have expected the existence 
of this solution, even before solving the equation for /i(t?), because the 
perturbation of this form can be generated by a coordinate transformation 
which does not violate our choice of synchronous coordinate systems and, 
hence, does not destroy our initial form (9) of the perturbed metric. (One 
can easily check that the function (32) is indeed a solution to Eq. (31).) 
Concretely, one can perform a small coordinate transformation 



a l drj 



In terms of new coordinates 77, x l the transformed components (9) take on 
the form 



g OQ = -a 2 (77) , 
where 



g o i = 0, g ik = a 2 (r]) 



h + C-, 
a 



(1 + hQ)S ik - hin 2 Q,i,k 
hi = h\ + Cn 2 j a~ l dr/ . 



(33) 
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The same transformation should be applied to the components of the 
energy-momentum tensor. Even if the original h, hi are zero, the trans- 
formed h, hi are not zero. After erasing the overbars in the transformed 
components of the metric, one returns to Eq. (9). The freedom of choosing 
different freely falling coordinate systems and corresponding spatial slices 
77 = const gets represented in the form of freedom to choose different solu- 
tions from the family of all solutions for the perturbations. All choices of C 
are equally well "physical" . The integral in Eq. (33) produces an additional 
integration constant which reflects the possibility to make a purely spatial 
transformation and to shift hi by a constant value, but we will not actually 
need this coordinate freedom. It follows from Eq. (33) that there are two 
functions (and many algebraic and differential combinations constructed 
from them) that do not contain C at all: 

u = h' + ajh , v = h'i n 2 h . 

a 

The solution (32) allows one to reduce the third-order differential equa- 
tion (31) to the second-order differential equation. To reach this goal we 
use the function it (77): 

u = ti + cry/i . (34) 

Obviously, this function vanishes on the solution (32). By substituting 
Eq. (34) into Eq. (31) we derive the equation for u(rj): 

u" + u'\2aj - ^ j +u 




n 2 - 2a' - a— - — = . (35) 



Note that the coefficients of this differential equation depend exclusively 
on the scale factor and its derivatives. This fact is a manifestation of the 
underlying interaction of the perturbations with the cosmological pump 
field. No special assumptions about the shape of the potential V(ip) or 
such things as "nonsimultaneous rolling of the scalar field down the hill" 
have been made whatsoever. 

Our next move is to transform Eq. (35) to the form similar to the equa- 
tion for gravitational waves. This will allow us to use certain results derived 
previously for gravitational waves and rotational perturbations [17,2]. In 
order to get rid of u' we introduce the function /j(ry) according to 

«\/7 , q «a 
u = — - — n . (36) 



It follows from Eq. (18) that the function 7(77) is nonnegative if the scale 
factor is governed by the scalar field (1) which we study here. However, 
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Eq. (35) is formally applicable to negative 7 as well. It may happen (as 
the author thinks) that Eq. (35) has a wider domain of validity and can be 
used, for other models of matter, with negative 7 too. If this is the case, 
one is free to modify Eq. (36) by using vTtI instead of ^7. Anyway, with 
the help of Eq. (36) one derives the equation 

A*" +V[n 2 -U(r))] 







where 



U (v) 




2 1 / 1 7 
a + a + a h 

7 



n" I 

a +a = — , C/i(r/J = -r 
a 7 Z 



C/ (r/) + C/i(r?) 



«77 



\i 2 + ^77" 



(37) 



The effective potential U(rj) can also be written as U(rf) = (a v / 7)"/ a v / 7 
which reduces our basic equation to the form 



n 2 - 







(38) 



(The function a^/7 can be related with the function z discussed in [18], see 
also the early papers [27].) 

Let us recall [3] that in the case of gravitational waves the potential 
U(rj) consists only of the Uo(r]) term, so that the basic equation is 



n 



a 
a 







(39) 



The gravitational wave potential Uq depends only on the first and second 
time derivatives of the logarithm of the scale factor: (In a)', (In a)". The 
potential U(rj) for density perturbations is more complicated and includes 
also (In H)', (In H)" , and (In H)'" . We note, however, that the potentials 
are exactly the same, and, therefore, the basic equations and solutions for 
density perturbations and gravitational waves are exactly the same, if 7 is 
constant, that is for the scale factors (6). For this class of pump fields, the 
general solution to Eq. (39) can be written in the form (for non-half-integer 



Kv) = (nr/) 1/2 



AiJpMnr)) + A 2 J_ iJ3+ i ) {nri) 



(40) 



Having two linearly independent solutions to Eq. (38) one can construct 
h(rj) and, hence, to find the rest of functions describing density perturba- 
tions. It follows from Eqs. (34) and (36) that 



ti 



7 . 

-jali H 1 — n 



(41) 
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and 

h( V ) = - [\^d V +-d (42) 

where r/o is some initial time where the initial conditions are to be imposed. 
The constant C entering Eq. (32) is denoted Cj at the i stage and will have 
the labels e and m at the e and m stages. All (complex) solutions to our 
perturbation problem for a given wave vector n are completely determined 
by three arbitrary and independent (complex) constants. Two of them de- 
fine a solution to Eq. (38) (these constants are A\, A2 when Eq. (40) is 
applicable). The third constant, Cj, describes the remaining freedom in our 
choice of coordinates. This remaining freedom is not a misfortune of the 
theory. On the contrary, it will later allow us to join our coordinate system 
right to the comoving synchronous coordinate system at the m stage. 

One can show by using Eqs. (24), (12), and (41) (and assuming ip' ^ 0) 
that 

r l 



<Pi(v) = -7= \fl h 
One can also find with the help of Eqs. (11), (13), and (41) that 



(43) 



Pi 
ei 



(44) 



// - fi \ 4a + j + 3(2^/7 a/i 



The quantity q, where cf /c 2 = pi/e±, plays the role of the velocity of sound 
for the high-frequency scalar field oscillations (see also Sec. VII). 

Similarly to what is true for gravitational waves, solutions to Eq. (38) are 
different for the high-frequency and low- frequency regimes. In the former 
case, n 2 3> |J7(r/)| and /x ~ e ±mr] . In the later case, n 2 <^ \U(r])\ and two 
independent solutions are /ii ~ a-y/j, H2 ~ ci^fy j drj / (a^/j) 2 . The functions 
^1, 1^2 generalize the corresponding solutions for gravitational waves by 
replacing a with Specifically for the scale factors (6), the solutions 

fii, /i2 are /ii ~ r] 1+ @ and /i2 ~ V~^i m agreement with Eq. (40). 

In the high-frequency regime, the term // dominates the other terms in 
Eq. (44). As one could expect, in this regime, the velocity of sound is equal 
to the velocity of light, cf = c 2 . In the low-frequency regime, that is when 
a given mode enters the under-barrier region, the dominant solution is 
(for a review, see [19]). Using this solution in Eq. (44) one can show that, in 
this regime, pi/e\ « q((3), that is the "velocity of sound" is the same as the 
one defined by Po/e' . In particular, pi/ei goes to —1 for the low-frequency 
scalar field solutions at the De Sitter stage, (3 = —2. 



17 



Equations (31), (35), and (38) have been derived under the condition 
(p ; / 0. However, the final formula (42) gives the correct result 
h{rf) = Ci/lo = const in the De Sitter limit 7 = 0, <p' = 0. One can 
analyze this case separately, referring to the starting Eqs. (23)- (26). One 
can see that Eqs. (24) and (12) give £' = 0, h = const. Equations (13), 
(25), and (23) give p\ = 0, V t<p = 0, e± = 0. Finally, Eq. (11) requires 
h[ = —7]n 2 h = —7]n 2 Ci/lo. But this solution for h, h[ is precisely solution 
(33) which can be eliminated by a coordinate transformation. In the De 
Sitter case governed by the scalar field (1) there is no density perturbations 
at all. Note that the function (pi(rj), Eq. (22), remains arbitrary and the 
wavelengths of these fluctuations are growing in the course of expansion. 
If one wishes, one can attach to ip\ such words as "inflation is pushing the 
waves beyond the De Sitter horizon". Nevertheless, the result will be zero, 
as long as (p± is not accompanied by perturbations of the gravitational 
field. This is an instructive example in order to realize that to "stretch 
the waves outside the causal horizon" is not all we need for generation of 
density perturbations (likewise, it is not sufficient to simply stretch the 
electromagnetic waves "beyond the horizon" in order to generate photons). 

In order to receive further insight into the statements discussed on pre- 
vious pages, the reader of these School Proceedings may be interested in 
performing some exercises. In particular, one will be able to evaluate cer- 
tain claims which exist in the inflationary literature and which gradually 
became "standard". 

It should be clear from the previous discussion that the essence of the 
parametric (superadiabatic) amplification mechanism is the interaction of 
the perturbations with the variable gravitational field of the FLRW uni- 
verses. It was emphasized long ago [3] that the condition of equality of 
the wavelength of the perturbation to the Hubble radius is, at most, a 
necessary condition for amplification, but not a sufficient condition. [This 
condition is not even a necessary condition if the velocity of sound of the 
density or rotational perturbations (in contrast to gravitational waves) is 
much smaller than the velocity of light. In this case, the amplification may 
take place long before the wavelength of the perturbation becomes compa- 
rable with the Hubble radius. What is a necessary condition for parametric 
amplification is that the frequency of the pump field should be compara- 
ble with the frequency of oscillations, not the wavelength of oscillations be 
comparable with the Hubble radius.] Even if the wavelength of a given 
mode has become comparable with the Hubble radius and the frequency of 
oscillations has become comparable with the frequency of the pump field, 
it does not mean that the amplification of the mode must take place. The 
wavelength of a given mode is always proportional to the scale factor a(ry), 
and the Hubble radius is always proportional to a 2 /a'. These are simply 
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their definitions and there is no indication of any interaction in these kine- 
matical facts. What we are really interested in is how the amplitude, not 
the wavelength, of the mode behaves. We need to know whether or not 
there is a coupling of the mode to the pump field, and what is the explicit 
form of the coupling (or, equivalently, what is the effective potential U (77) 
in the dynamical equations such as Eq. (38) or Eq.(39)). The wavelength 
and the Hubble radius may "cross" each other in course of their evolution 
as many times as one wants, it does not guarantee any superadiabatic am- 
plification, if there is no appropriate coupling. However, in the inflationary 
literature the reason for amplification is usually seen in that the inflation 
"stretches" the waves and acts as a "magnifying lense" or as a "micro- 
scope". The explanation of the generating mechanism goes usually along 
the lines: "oscillations inside the horizon" , "stretching" , "crossing outside 
the horizon", "freezing", "crossing back" and ... the created perturbations 
seem to be ready for comparison with observations. In order to appreci- 
ate the difference between the kinematics and dynamics one can make the 
following: 

Exercise 1. 

Take your favorite inflationary model, whichever you like. Consider two 
test waves with the same wavelength: one is gravitational wave, another 
is electromagnetic wave. Show that their wavelengths behave in exactly the 
same manner ("crossing, stretching, crossing back"). Show that the gravi- 
tational wave will be superadiabatically amplified (and, hence, can be gen- 
erated quantum-mechanically) but the electromagnetic wave will not (the 
corresponding potential in the dynamical equation analogous to Eq. (39) is 
identically zero). 

Let us now turn specifically to density perturbations in models governed 
by the scalar field (1). The archetype inflationary model is the exponential, 
or De Sitter, expansion which corresponds to the case (p' a = 0, 7 = 0. This 
is the model which generates the scale-invariant spectrum of cosmological 
perturbations. Here and in the rest of the paper we will be interested in 
what theory says about the numerical values of the amplitudes that form 
the scale-invariant (or almost scale-invariant) spectrum of density pertur- 
bations. It was argued in the text above that, despite the presence of the 
nonvanishing scalar field perturbations, there is no density perturbations 
at the De Sitter stage. The scalar field perturbations may oscillate and 
then "cross the horizon" but there will be no accompanying density pertur- 
bations neither "inside" nor "outside" the Hubble radius. In other words, 
in the De Sitter case, the solution to the linearized Einstein equations is 
a purely coordinate solution totally removable by an appropriate coordi- 
nate transformation. For models that slightly deviate from the exponential 
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expansion (the quasi-De Sitter expansion) the function 7(77) is no longer 
strictly zero and some amount of density perturbations is present. How- 
ever, if 7 is a small number the amplitude of density perturbations is also 
small, according to Eq. (42), regardless of whether the wavelength is shorter 
or longer than the Hubble radius. Now, a considerable part of papers on 
the subject of perturbations deals with the so-called gauge-invariant poten- 
tials [16]. In order to make contact with other studies one can perform: 

Exercise 2. 

Show that the gauge-invariant potentials for density perturbations are 
strictly zero in the De Sitter case (7 = 0) and they are small if 7 is a small 
number, 7 <C 1. 

The next Exercise will be especially instructive in view of the previous 
ones. I believe that the lack of clear physical picture of the phenomenon 
may lead to the loss of immunity to wrong conclusions. It follows from 
Exercise 2 that the amplitude of density perturbations was very small if 
the expansion rate was very close to the De Sitter expansion at the time 
when the wavelength of our interest was first "crossing outside" the Hubble 
radius. Let the wavelength of our interest be such that it "crosses inside" 
the Hubble radius today. What is the amplitude of density perturbations 
in this mode today? The "standard" inflationary formula says that 



Sp 
P 



H 2 



H <K*1 



where (j) is supposed to be taken at the time (t\) when the wavelength of the 
perturbation "crossed outside" the Hubble radius during the De Sitter or 
quasi-De Sitter epoch (see also Introduction to this paper). In other words, 
according to this formula, the predicted amplitudes of the scale-invariant 
spectrum go to infinity. 

There exist many "intuitive" derivations of the "standard" formula. 
They can only illustrate the fact that the physical intuition of different 
people is different. Apparently, it is perfectly in accord with the intuition 
of some theorists that the amplitude which was almost zero initially has 
transformed by today into almost infinity, simply because the wavelength 
of the perturbation was first crossing the Hubble radius when the expansion 
rate was almost the De Sitter one. Instead of questioning the validity of this 
formula, theorists went on deriving the far-reaching consequences of it. This 
formula is essentially the basis for discarding certain scalar field potentials 
on the grounds that they generate "too much" of density perturbations and 
for claims that the contribution of gravitational waves to the quadrupole 
anisotropy is "negligibly small" in the limit of the De Sitter expansion. 
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For instance, the authors of the recent paper [28] say that "In the limit 
that n s — > 1 (p — > oo) gravitational waves do not contribute to COBE's 
signal, and the power spectrum for zeta is the flat, Zel'dovich spectrum". 
Obviously, there is nothing wrong with the generation of gravitational waves 
in the De Sitter limit (p — ► oo). The reason for the above conclusion is the 
authors' formula (5.15) (the analogue of the "standard" formula) which 
states that the amplitudes of density perturbations today are arbitrarily 
large in the De Sitter limit (the denominator goes to zero). One can hear 
sometimes that it would be a disaster if the gravitational wave contribution 
to the large-angular-scale anisotropy was dominant since there would be no 
scalar perturbations to seed galaxies. I, personally, think that it would be 
a disaster to build serious conclusions on unreliable theoretical formulas. 

Since the intuition of different people is different, we will try to see how 
the "standard" formula was derived from equations. We will discuss two 
papers [18,29]. [I don't want to be unfair with respect to these authors: we 
will discuss these papers precisely because they are discussable, they contain 
the formulation of the problem, equations and certain logic in calculations. 
In my opinion, this cannot be said about many other papers on this subject.] 

The basic equation of Bardeen et al. [29] is their equation (2.23), the ba- 
sic equation of Mukhanov et al. [18] is their equation (6.49). These equations 
were obtained by manipulating with the original perturbed Einstein's equa- 
tions in order to express them in terms of the gauge-invariant potentials. 
Both groups of authors introduce the gauge-invariant quantity zeta and use 
it for derivation of the final conclusions. The quantity zeta by Mukhanov et 
ol-, Cm.f.b, is defined by their equation (6.66), the quantity ( by Bardeen et 
a l-> Cb.s.t.5 is defined by their equation (2.45). Both group of authors neglect 
in their equations the term which contains the square of the wave number 
and work with the approximate equations. They claim that in this approx- 
imation the quantity £ is constant. Bardeen et al. say that Cb.s.t. "is nearly 
time independent throughout the period when k/SH <C 1". Mukhanov et 
al. say that Cm.f.b is conserved "when considering wavelengths far out- 
side the Hubble radius for which V 2 <1> can be neglected". In particular, 
Mukhanov et al, by using their conservation law, arrive at their equation 
(6.67) which says that that the final amplitude of perturbations can be 
arbitrarily large if the initial equation of state (at the first Hubble radius 
crossing) was sufficiently close to the De Sitter case. [There is a misprint 
in their formula (6.67): the position of quantities w(ti) and w(tf) should 
be interchanged.] On the basis of their formula (6.67) Mukhanov et al. 
come to formula (6.68) which is an analogue of the "standard" formula 
and which is "in agreement with (6.66)" and "with the results of previous 
investigations ..." [They give references to the papers listed as [29,30] in 
the present paper.] 
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Exercise 3. 

Using the definitions of the gauge-invariant potentials and notations of 
the present paper show that the basic equations of Bardeen et al., their 
(2.23), and of Mukhanov et al, their (6.49), are totally equivalent to the 
equation 
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Show that this equation can be reduced to the form 
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a 2 7 
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(38') 



where X is an arbitrary constant. Thus, the basic equations of Bardeen et 
al. and Mukhanov et al. are totally equivalent to equation (38'). Now, using 
the definitions of the zeta functions show that 



,M.F.B. 



1 1 

2n 2 a 2 7 



a 2 7 



and 



Cb.s.t. 



1 

2^2 



a 2 7 



a 2 7 



+ 



6a \ a^/7 



so that Cm.f.b. and Cb.s.t. coincide in the limit n 2 — ► 0. Show that the 
approximate equation which Bardeen et al. and Mukhanov et al. work with 
corresponds to the omitting of the term with n 2 in equation (38'). Show 
that equation (38') gives in this approximation 



Cm.f.b. 



X 
2^2 



which is a constant indeed. And, finally, compare the exact equation (38') 
with the exact Eq. (38) of the present paper to find out that the constant 
X must be equal to zero. Thus, in the considered approximation, the con- 
servation law C(i?) = C(*/) translates into the statement = 0, and the 
claimed "standard" formula does not follow from it. 



One additional point: note that the condition of applicability of the 
approximate equations (without the term with the square of the wave num- 
ber) is 

(aV7)" , 



n 2 < 



a\/l 
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(loosely speaking, the wave is under the potential barrier), not n 2 <C a 2 
(the wave is outside the Hubble radius). This will be important in the next 
Section. 



4. Late Time Evolution of the Perturbations at the Initial Stage 

The main uncertaintes about the evolution of the very early Universe refer 
to the times that preceded the epoch of the primordial nucleosynthesis. 
Whatever was the initial stage, even if it was governed by the macroscopic 
energy-momentum tensor of the collection of superstrings, it supposedly 
went over by that epoch into the radiation dominated stage governed by the 
scale factor a(rf) = l a e (r) — r) e ). The constants a e , r) e are to be determined 
from the continuous joining of a{rj) and a!{rj) at the time r\ = r)\ of transition 
from the i stage to the e stage. If the % stage is described by the scale factors 
(6), one derives 

a e = -(1 + (3)\m\ 13 , Ve = Y^TpVi ■ 

In further applications, we intend to use simple solutions (40), (42) and 
to make their appropriate joining with perturbations at the e stage. How- 
ever, the function 7(77), being equal to the constant 7 = (2+/3)/(l+/3) at the 
i stage, and to the constant 7 = 2 at the e stage, experiences a finite jump 
at the transition point 77 = 771 . This presented no problem for gravitational 
waves, since 7' (7?) did not enter the gravitational wave potential Uq{t]). But 
this becomes important for density perturbations, since the U\(rj) part of 
the potential acquires increasingly growing values at the end of the % stage 
for steeper and steeper transitions of 7(77). 

To deal with the problem, we introduce a parameterized set of smooth 
functions 7(7?) that approximate the step function in the limit of the pa- 
rameter e going to infinity: 

4 _|_ 3/J R 

^ ) = WW) + WW) ianh[e{ll ~ 11l)] ■ (45) 

For large negative values of 77 the function (45) goes to (2 + (3)/(l + /?), and 
for large positive values of 77 it goes to 2. We may surround the transition 
time 77 = r]i by a thin "sandwich" with boundaries at 771 — a and 771 + a. 
The asymptotic values of 7(77) are already reached with arbitrary accuracy 
at the boundaries, if e is sufficiently large, e 3> l/cr. 

The function (45) can be integrated, see Eq. (4), to produce the function 
0(77): 



a(rj) 1 + 2(1 + 13) 



I e e(n-m) + e -e(v-vi) 
77 - 771 + - In 



(46) 
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Again, for sufficiently large e, the function a(rj) quickly approximates 
(1 + P)/rf to the left of the transition point, and 1/(77 ~~ We) to the right of 
the transition point. These are the values of a(rj) that are appropriate for 
the i stage (6) and the e stage, respectively. 

The divergent functions 7', 7 /2 and 7" participate in the potential 67(77) 
(Eq. (37)). The function 7' grows as e at the point 77 = 771. The function 
7" is equal to zero at 77 = 771, but it grows as e 2 slightly to the left of this 
point, and it grows as — e 2 slightly to the right of this point. We assume 
that the transition to the e stage has completed at 77 = 771 + a, and we 
let a to go to zero. The function U\(rj) compresses and stretches to the 
arbitrarily large positive and negative values when a goes to zero and e 
goes to infinity. In particular, the potential U(rj) goes through zero and the 
condition n 2 <C |(ay / 7)"/( a \/7)l can be violated even if the wavelength is 
still longer than the Hubble radius. Examining Eq. (38), one can expect 
that the value of /i'(r/) at 77 = 771 +0 will be different from the value of //(r/) 
at 77 = 771 — a. The integration of //" in the limits from 771 — a to 771 + a gives 
a jump in fj! which depends on the value of the integral from the divergent 
part of the potential U\(jj). Fortunately, it is not fJ{w) itself, but a particular 
combination (■ v /7/a)[/t / — M( a + t'/^t)] that we will need to know in our 
further calculations. This simplifies the analysis. Due to Eqs. (29) and (41) 
this combination is precisely the function 77(77) introduced in Sec. III. 

In terms of the function v(rj), where 



v = - — 
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" ] , (47) 
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the basic equation (38) takes on the form 

(a 2 v)' = -n 2 a^j n . (48) 

The integration of this equation over the thin "sandwich" shows that 
v\ v=m +o = v\ v=m -o- In other words, the function // — tj(q+|^-) taken right 
at the beginning of the e stage is equal to the value of this function taken 
right at the end of the i stage (6) times the factor ^ \Jl+@- ^ n addition to 
the conditions: 



7l?7=?n-0 — y j Q , 7l??=?n+0 — ^ , A*|tj=tji— — H??=r/i+0 j 

this establishes the rules for going through the "sandwich" . 

The meaning of this procedure is the following. Instead of following the 
gradual growth of the function 7(77) from some initial small value up to the 
final value 7 = 2 (in which case it is difficult to find an exact solution to 
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Eq. (38)) we concentrate the entire evolution of 7(7/) to a short interval of 
time. In the limit, we consider a sharp jump of 7 from one constant value 
(which may be zero as at the De Sitter initial stage) to another constant 
value, but we take into account the form of the potential U{rj) and the 
change in the initial exact solution which is generated by this jump. 



5. Density Perturbations in the Perfect Fluid Matter 



We will now consider Eqs. (11)-(14) at the perfect fluid stages governed 
by the energy-momentum tensor (15). Similarly to the scalar field case, 
the longitudinal-longitudinal part of stresses vanishes, p\ = 0. For the eas- 
ier handling of arbitrary e\, p\ one can introduce the following notations: 
Pi/ei = cf /c 2 and p' /e' = c 2 s /c 2 , see Eq. (7). These definitions are conve- 
nient but, generally speaking, they have only formal meaning, since both 
Pi/e\ and Pi/ e o can De negative. However, in certain regimes, the quantity 
c\ is a genuine longitudinal velocity of sound (see Sec. VII). Our first in- 
tention is to derive the equation for h(rj), analogous to Eq. (31) and valid 
for arbitrary nonzero pi/e±. Specific cases cf = |c 2 and cf = will be 
considered separately. 

In order to derive the equation for h(rj) one can essentially repeat the 
steps that have lead to Eq. (31). Find h\{r]) from Eq. (11) and plug it 
into Eq. (14). Use Eq. (13), the first derivative of this equation, and the 
definition of cf/c 2 . As a result, one arrives at the equation 
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Now, introduce the function u{rj) according to Eq. (34) and use it in 
Eq. (49). Equation (49) can be reduced to 
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(50) 



The last term in this equation vanishes if 
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which integrates to 



c 2 const 
cf aa 2 7 



An assumption which is usually made for perfect fluids is: 

cf = c 2 s . (52) 

Note that Eq. (52) is certainly true for matter with the equation of state 
p = qe, where q is a constant, but Eq. (52) is not true in general, and it 
is not true for the scalar field matter (1) (unless one considers the under- 
barrier region where Eq. (52) is true approximately, see Sec. III). Due to 
Eq. (51) the last term in Eq. (50) cancels out. (If we have not assumed (51), 
the function h(rj) = C(a/a) would not have been a solution to Eq. (49).) 

We can now introduce the function v(rj) according to (compare with 
Eq. (36)): 

u = ^ c s v . (53) 
a 

In terms of v{rj), Eq. (50) takes on the form 

= , (54) 
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The potential W(rj) depends exclusively on the scale factor and its deriva- 
tives. Having a solution 1/(77) to this equation and using Eqs. (53) and (34), 
one can construct the function h{rf) and the rest of perturbations. Similarly 
to the scalar field case, all solutions for perturbations with a given n arc 
defined by three constants one of which describes the remaining coordinate 
freedom. These constants are expressible in terms of the constants given 
at the preceding i stage through the joining of the perturbations at the 
transition time 77 = 771 from the % stage to the perfect fluid stage. 

We will now consider Eq. (54) specifically at the radiation-dominated 
stage p = |e. We have cf = c 2 = (l/3)c 2 , 7 = 2, and the scale factor 

a(rj) = l a e (rj - rj e ) (55) 

where the constants a e , rj e are to be determined from the continuous joining 
of a(rj) and a' (77) at the transition time rj = r)\. The potential W(rj) vanishes. 
Equation (54) simplifies to the familiar equation 

v" + -n 2 u = (56) 
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which describes the time-dependent part of sound wave oscillations in the 
radiation-dominated fluid. In what follows, we will be using the general 
solution to this equation written in the form 



Bie -^-^) + B2e ^(v-v.) 



(57) 



where B\, B2 are arbitrary and independent (complex) numbers for each 
individual wave vector n. 

The function h(rj) is determined by a known solution for v{rj) and a 
coordinate solution with arbitrary constant C e : 



h(v) = ~ u dr) -\ — C e 
a J rn a 



All other functions are expressible in terms of h(rj). In particular, 

h\ = -[Sh" + 9ati + n 2 h] . 



(58) 



(59) 



The general solution (57) is always oscillatory in ry-time. The perturbed 
energy density, pressure, and the associated gravitational field h(rj), hi(rj) 
oscillate in space and time as they should do for sound waves. However, if 
one considers these oscillations at intervals of time shorter than their period, 
they may appear as consisting of "growing" and "decaying" solutions. In 
particular, this happens if one considers the relatively long waves, 

~ We) = 2/1 *C 1, at their first oscillation since the beginning of the 
e-stage, that is while the condition -^(w ~ We) "C 1 is satisfied. Under this 
condition, the function h(rj), Eq. (58), can be approximated as 



h( V ) = ^ + ^ + B 2 + 

7] z 7] 



(60) 



where 



and 



^(V-We), 



c e + 



Si = 



B1 + B2 



Bo = 



-vnjBi - B 2 ) 
2\/3l Q a e 



iV3 



n 



Bi 1 



-iyi 



Bo 1 



The common practice [15] is to use the coordinate freedom for elimi- 
nation of the "most divergent" term in Eq. (60), which is also the "most 
decaying" term, if one goes forward in time. This is achieved by such a 
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choice of C e that C e = 0, and the first term in Eq. (60) vanishes. Then, the 
energy density perturbation 5e/eo (use the definition 

Se/e = ^Lq 

and calculate ne\ according to Eq. (11)) can be approximated as 

5 S = ^B 1 f)+ 1 -n 2 B 2 n 2 + --^Q . (61) 

The part of Eqs. (60) and (61) which depends on the coefficient B\ is 
usually called the "decaying" solution, while the part with the coefficient 
B 2 is called the "growing" solution. 

Despite the possibility of identifying (quite artificially) the "growing" 
and "decaying" solutions, density perturbations at the e-stage form a col- 
lection of traveling sound waves with arbitrary amplitudes and arbitrary 
phases, as long as constants B\, B 2 are arbitrary and independent. One 
should not think that simply because the sound waves have spent some 
time "beyond the horizon", nf « 1, they would transform into standing 
waves at later times of their history after they came "inside the horizon" , 
nfj 2> 1. To illustrate this point, let us take into account the spatial part of 
the perturbations and consider the contribution h n (rj,x) of a given mode 
n to the total field h(rj, x) = J2n h n (n, x). This contribution can be written 
as 

fcofox) ~ h n e inx + h* n e- in * 

% n e~ inf] - B 2n e inf} ) e inx + (sj n e inf > - B* 2n e~ in ^ 
= 2|Bi n | cos(n?7 — nx — ipi n ) — 2\B 2n \ cos(n^ + nx + </?2n) (62) 

where B ln = l^inle^ 1 ", B 2n = |B 2 n|e i,Wn . The last line in Eq. (62) shows 
explicitely that, in general, one is dealing with waves traveling in opposite 
directions with arbitrary amplitudes and arbitrary phases. A standing wave 
can only be "generated by hand", by assuming that the constants B\ n , 
B 

2n are strictly related. This happens if one declares that he/she is only 
interested in the "growing" solution and puts B\ = 0. Then, the complex 
amplitudes Bi n , B 2n become related: |i?2n| = (— l) fe+1 |-Bin|j 
9?2n = fin — kir, (k = 0,1,2,...), and the last line in Eq. (62) can be 
transformed to 

h n (r], x) ~ 4|5i n | cos nfj cos(nx + tpi n ) 



which is a standing wave indeed. 
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Standing sound waves at the e stage are responsible for so-called 
Sakharov's oscillations [20] in the power spectrum of density perturbations 
in the present Universe, at the m stage. As we have shown, standing sound 
waves cannot originate somehow automatically at the e stage, simply be- 
cause of the transition from the "growing" / "decaying" regime to the oscil- 
lating regime (see also [21]). If one works with classical density perturba- 
tions at the e stage and makes no additional assumptions, one can say noth- 
ing about the necessity of standing waves, except of postulating this. The 
point of this discussion is that the quantum-mechanical generating mech- 
anism, which we are considering in this paper, does really create standing 
waves. Standing waves arise for gravitational waves, rotational perturba- 
tions, and density perturbations. The physical reason for this is that the 
waves (particles) are generated in correlated pairs with equal and oppo- 
sitely directed momenta (the two- mode squeezed vacuum quantum states). 
This is true for waves of any wavelength, as soon as conditions for their 
generation are satisfied. Technically, as we will see later, the second term in 
Eq. (60) taken at the beginning of the e stage turns out to be much smaller 
than the third term in Eq. (60), that is B± + B2 ~ for y\ <C 1. 

We should now discuss a great difference between sound waves and 
gravitational waves with regard to their evolution in time. The velocity of 
sound waves at the e stage is only y/3 times smaller than the velocity of 
gravitational waves, but their amplitudes behave drastically different. The 
amplitude of a gravitational wave decays as a' 1 in course of time. The 
amplitude of a sound wave, as one can see from Eq. (58), decays as a~ 2 , 
since a ~ a" 1 . This leads to a difference in solutions even for relatively 
long waves which did not complete even one cycle of oscillations during 
the entire e stage from rj = rji to r/ = r/2- As an illustration, let us consider 
sound waves which barely reached the oscillating regime by the end of the e 
stage. Their wave numbers satisfy the condition -%(f?2 — Ve) = zr = V2 ~ 1- 
These are the waves whose wavelength was of the order of the Hubble 
radius at the time of transition from the e stage to the matter dominated 
m stage. Assuming that the present day Hubble radius is Ih ~ 6 • 10 3 Mpc 
and that the present day scale factor a(r//j) is cl{t]r) ~ 10 4 (1(7/2)) their 
wavelength today A c = 2 r na(j]u)ln c is about 220 Mpc. The usual practice, 
in addition to eliminating C e , is to concentrate on the "growing" solution, 
that is to assume that at the beginning of the e stage the second term in 
Eq. (60) is smaller, or at least not larger, than the third term. Under these 
conditions, the final numerical value of the function ^(77) is of the same 
order of magnitude as the initial value, hfa) ~ h(rji), for the wavelengths 
of our interest. In other words, if the preceding i stage has produced h(rj) 
with some initial numerical value h(r)i), this number will effectively be 
transmitted to the beginning of the m stage. 
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However, in the very same coordinate system where the "most decaying" 
term in Eq. (60) was eliminated, the time derivative h'(r)) was left large. 
The final value of h'(rj) for the waves of our interest is h!(r]2) ~ n/i (772). 
According to Eq. (12), the function h'(rj) describes velocity of matter. This 
velocity will be inheritted by matter at the matter-dominated stage. But 
this is not velocity of the fluid elements with respect to each other, this 
is not velocity describing deformations of the medium, and the functions 
h'(rj), h'^r]) are not the ones that we may use for our later calculation of the 
varations in CMBR. The function /i'(r?) describes velocity of matter with 
respect to the coordinate system that we have chosen for our convenience 
of eliminating the "most decaying" term. At the m stage, however, we are 
more interested in the comoving coordinate system. We are interested in 
the density contrasts and deformations of the fluid itself, we are interested 
in the components of the accompanying gravitational field that we may 
use for calculations of 5T/T. We need a coordinate system which provides 
vanishing of h'(r]) by the end of the e stage. This requires a different choice 
of C e . Under this new choice of C e , the first term in Eq. (60) survives, 
and numerically the same, as in the previous example, initial value h(r]i) 
transforms into a small number h(r]2) ~ yf/i(r/i) by the beginning of the 
m stage. This is what we need to keep in mind when we will compare the 
amplitudes of density perturbations and gravitational waves. 

Finally, let us consider density perturbations at the m stage, p = 0. At 
this stage, one has cf = c 2 = 0, 7 = |, and the scale factor 

a(rj) = l a m {v ~ Vm) 2 • (63) 

The scale factor and its first time-derivative are continuous at the time 
r] = r/2 of transition from the e stage to the m stage. Therefore, 
o m = a e /A(r]2 — i] e ), f]m = — f]2 + 2r/ e . It follows from Eqs. (13) and (14) 
that the general solution for h(r)), hi(rj) has the following form 

h = C\ H — C m , b! = — — — C m , 
a 2 a 

h\ = \cm\ri - V m) + -n 2 C m + C 2 ^ ~ Vm) " . (64) 

The energy-density and velocity perturbations can be found from Eqs. (11) 
and (12). Similarly to the preceding i and e stages, the perturbations are 
completely determined by three constants C\, C2, C m one of which, C m , 
reflects the remaining coordinate freedom. 

The constant C m is entirely responsible for a possible relative velocity 
of our fluid with respect to a chosen synchronous coordinate system, see 
Eqs. (12) and (10). A given coordinate system is not comoving, T % Q / 0, as 
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long as C m 7^ 0, h' ^ 0. But we know that for a dust-like fluid without 
rotation one can always introduce a coordinate system which is both syn- 
chronous and comoving. This is reflected in our ability to remove the C m 
term from h(rj) and hi{r]) by a coordinate transformation (33). Thus, the 
choice C m = in Eq. (64) is not a restriction of the physical content of the 
problem, it is an allowed choice of the coordinate system. 

It is here that we will eventually restrict our coordinate freedom, we 
will put 

C m = . (65) 

The constants C e and Q will not be arbitrary any longer, they will be 
determined from the continuos joining of solutions. We need the comov- 
ing coordinate system for simple and appropriate formulation of the 5T/T 
problem. We are interested in the temperature of CMBR and its anisotropy 
seen by a comoving observer, that is by an observer whose world line is one 
of the matter's world lines. One of these idealized comoving observers is 
an observer on Earth (up to accuracy of some nonzero peculiar velocity, 
local rotation, etc.). We are much less interested in feelings of an observer 
who wanders in the Universe with arbitrary time-dependent velocity. In 
the comoving coordinate system, the world line of a comoving observer is 
described by simple equations x l = const. 

Upon the choice of C m = 0, the perturbations reduce to 



\3 



M,) = - c in ^-^--c 2 ^-^ . (60) 

From Eqs. (66) and (11) one can derive the familiar expression [15] 



Se 



1 n 2/ \2 l n (V2-Vm) 3 



(67) 



(One may wish to correct a misprint in Eq. (115.21) of Ref. [15]: the de- 
caying solution behaves as r]~ 3 , not as the printed r/~ 2 .) As long as the 
constants Ci, C2 are arbitrary, the power spectrum of the density perturba- 
tions is arbitrary. In particular, there is no Sakharov's oscillations, a priori, 
and they do not arise simply because of the transition from the e stage to 
the m stage. For instance, one can start from a perfectly smooth spectrum 
at 7] = 7j2 and extrapolate these data back in time up to the beginning of 
the e stage. 

For the further calculations of 5T/T we will need the first time deriva- 
tives of the gravitational field perturbations at the m stage in the comoving 
coordinates. Since for the scalar component of the perturbations one has 
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h' = 0, it is only the longitudinal-longitudinal component h\ that is effec- 
tive. 

6. Joining the Perturbations at the Three Stages 

We are now in the position to start our operation of joining the solutions at 
i, e, and m stages. We want to derive from the first principles the expected 
density perturbations at the m stage. Of course, the result will depend on 
the unknown behaviour of a{rj) at the i stage. But this is precisely why 
we are doing this study: we try to learn something about the evolution of 
the very early Universe by deriving the expected variations in CMBR and 
comparing them with the observations. 

The general rule for joining solutions to Einstein's equations is to match 
from the both sides the intrinsic and extrinsic curvatures of the transition 
hypersurface [22]. For our solutions written in the class of synchronous 
coordinate systems, this translates into the continuity of the spatial metric 
and its first time derivative. Since we have already assumed that a(rj) and 
a'(rj) join continuously, it is the continuous joining of h(r]), hi(rj), h'{rj), 
and h'^rj) that should be ensured. In fact, it is sufficient to follow h(rj), 
h! (77) , and h\ (77) as hi (77) is derivable from h\ (77) at all three stages up to the 
integration constant which can be removed by the remaining integration 
constant in Eq. (33) anyway. It is convenient to write h\ at the i and e 
stages, respectively, in the form (use Eqs. (29), (59) and (42), (58)): 



We will denote a(rj), a{rj) at 77 = 771 and 77 = 772 by cti, ol\ and 02, «2 
respectively. It is also convenient to introduce the parameter y = -^=(77 — 7? e ) 
and its values y±, yi at the transition points. 

We will first make the joining of solutions in general form, without 
adopting any particular coordinate system, any particular behavior at the 
i stage, and any particular wavelength of the perturbations. 

Let us start from the i-e transition. Whatever was the i stage and its 
late time behavior, it produced certain 7/(771), 7/(771) and ended at r] = 7/1 
with some values of the scale factor and its derivatives. For generality, we 
do not assume, for the time being, that 7(7/1) is exactly 2 and 7' (771) is 
exactly zero. 

From the continuous joining of h(rj), h'(rj), and ^(77) one derives: 




— h + -{u' - av) 



(68) 



a a 




(69) 
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B 1 e- iyi + B 2 e iyi = ^7 /J + a(2 - 7) 



H^dr] + Ci 



no 



B 1 e~ 1 ^ (1 + i m ) + £? 2 e iyi (1 - iyi) = - -a ^7 



/x - jU a + 
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(70) 



• (71) 



All functions in these equations are taken at 77 = 771 so that, for instance, 
7 means 7(7/1), // means /j'(r/i), 7' means 7' (771), etc. For a more compact 
record we will also use the following notations: 



ai v ai 



m'(?7i) - M(m 



1 7 ' 
«i + --(r/i) 



and 



— Cj + 



Equations (69)-(71) allow us to express the constants B±, B2, and C e de- 
scribing the perturbations at the e stage entirely in terms of the output 
values of the functions defined at the i stage. 

Let us now turn to the e-m transition. Again, from the joining of h{rf), 
h'(rj), and h\{rj) one derives 



3 



V3 



7? 



B x e~ iyi s + B 2 e iyi s* , 



(72) 



a 2 



Ci = --^C e + 



3a2 



6a 2 7/ 2 



[Si e~ tyi (s + 3t/2 e~ ld ) + B 2 e tyi (s* + 3y 2 e td )] , (73) 



C 2 



^^C e + ^{ Bl e-** [e-"(-5 + 27/1 " 6i W ) + 
5a2 5a2 1 



+ B 2e ^[ e ^(-5 + 27/ 2 2 + 6f7/2)-i7/ 2 ]}, (74) 

where d = y 2 — yi = ~^(V2 ~ Vi), s = 2i + (7/2 — 2i)e~ ld . Everything at the 
777, stage is known as soon as the coefficients C±, C 2 , and C m are known. 
They are expressed in terms of the coefficients B\, B 2 , and C e attributed to 
the e stage. Since these numbers, in turn, are known implicitly in terms of 
the coefficients attributed to the i stage, we have linked the very beginning 
with the very end. 

Our next step is to impose the requirement (65) and to choose the 
comoving coordinate system at the m stage. From Eqs. (72) and (69) one 
can find 

CiD = 2^^ai{a[3«i(l - iyi) + vi] - s*[3ui(l + i yi ) + «i]} , (75) 
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where 

D = 2yf + i?-^[s*{l + iy 1 )-s(l-iy 1 )} 

= 2yj + (2- 1 )[2-(2 + y 1 y 2 )cosd-{y 2 -2y 1 )smd]. (76) 

Equation (75) says how to choose Cj at the i stage in order to match right 
to the comoving coordinate system at the m stage. We can now put Eq. (75) 
into Eq. (70) and solve Eqs. (70) and (71): 

B ^ e ~ m {6yi(l-iyi)«i-[(2-7)**-2j/i]t;i}, (77) 

B 2 e iW = -7^— {%i(l + %)«!- [(2 - 7)« - 2yi]«i} . (78) 

Substituting these formulae and Eq. (75) into Eqs. (73) and (74), we reach 
our goal — the finding of C\, and C 2 in the comoving coordinates: 

d = - Gl {3ttij/i(sin d + yi cos d) - vi[(2 - j)(cos d - 1) - y 1 sin d}}, 
6a 2 OL\D 



(79) 



2oi 

C 2 = - ^^^ {3^1^/1 [(10 - 3yf - lOyu/2) sin d 



+ (-10y 2 + 10yi - 3yiy|) cos d] 

_ ^{-2(2 - 7 ) (5 + y\) - [yi(10 - 3y|) - 10y 2 (2 - 7)] sin d 
+ [(2- 7 )(10-3y 2 2 ) + 10y iy2 ]cosd}}. (80) 

So far, no approximations have been made. We will start making them now. 

The numerical value of the denominator D, and hence the absolute 
values of C\ and C 2 , depend critically on how close the exiting value 
7(771) =71 is to 2. We are interested in wavelengths that were longer than 
the Hubble radius at 77 = 7/1. Their wave numbers satisfy the requirement 
y± <C 1. On the other hand, y 2 /yi = a 2 /ai 3> 1, and d = y 2 — y\ ~ y 2 - 
The wavelenghts longer than the present day A c = 2ira(riJi)/n c w 220 Mpc 
correspond to small y 2 , and n < n c . These are the wavelenghts of the major 
interest for the discussion of the large-angular-scale anisotropy in CMBR. 
In the approximation of small y 2 , two leading terms in D are 

n ~ 2 , 2 - 71 4 
D w 7i yi + — 77 2 . 

The second term is much larger than the first one (and, hence, the expected 
Ci, and C 2 are hopelessly small) for all 



ai 
n c — 

a.2 



^/l2 7l /(2 -71) < n < n c 



34 



unless the exiting value 71 is so close to 2 that the second term can be 
neglected. In order to deal with the most favorable situation and not pro- 
liferate complications, we will assume that 71 = 2 and (7 / /7)( 7 ?i) = 0- We 
know, see Sec. IV, that the transition from the very end of the i stage (after 
a thin "sandwich" interval) to the very beginning of the e stage can be made 
arbitrarily smooth (at least, in theory). So, we will be using D « 2y\ <C 1. 

We should now take into account the fact that v\ <C u\ for all wave- 
lengths of our interest. Indeed, we are interested in modes that have inter- 
acted with the potential barrier in Eq. (38) and have been amplified at the i 
stage. For the scale factors (6), their wave numbers satisfy the requirement 
(nr/1) 2 <C (3{(5 + 1) which translates into the condition y\ <C /3/3(/3 + 1), 
or simply y\ <C 1. We will derive the approximate formulae valid in the 
leading order by the parameter y\. 

As we know, there are two independent solutions in the under-barrier 
region. Which of them dominates is the matter of choice of the initial con- 
ditions at r] = rj (choice of constants A\, and A2 in Eq. (40)). For classical 
solutions, one can choose the initial data in such a way that there will be 
no amplification at all, or there will be even attenuation instead of amplifi- 
cation. However, a "typical" choice of initial data at r\ = rj Q , which amounts 
to the averaging over the initial phase (or a rigorous quantum-mechanical 
treatment), always leads to the dominant solution jjl ~ a^/7, and to am- 
plification [19]. This is the choice that we imply here and will justify later, 
Eq. (102) in Sec. VIII. Since v ~ {^/a^f) 1 and jj, ~ a^/7, the quantity v\ 
is relatively small. Concretely, for solutions (40), one has 

(81) 

and 

T) A 1 

2/^+2 r(/3 + §) 

so that 

Ul 2/3 + 3 yi 

Thus, we can neglect all the terms containing v\ in Eqs. (77)- (80). 
In the leading order, one has 

Bi &-B 2 n^- V2n( m ) = B 
2yi 

B 1 + B 2 ^Byf . 
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These formulae ensure the standing wave pattern for all wavelengths at the 
e stage and the power spectrum modulation at the m stage (compare with 
Sec. V). The leading order expressions for C\, C2 are as follows: 

d » y/2 n(rn)wa. d (82) 

2a2 yi 



C 2 » -■= V2 n(rn) 

5a 2 yi 



(10 — 3yf) sin d — lCh/2 cos d 



(83) 



We will give a brief analysis of Eqs. (82) and (83). The growing and 
decaying components of hi(rf) and 5e/eo, see Eqs. (66) and (67), are of the 
same order of magnitude at rj = 772 for all wavelengths. The coefficients Ci, 
C2 are smooth for long waves, d ps 2/2 "C 1, n <C n c : 



2ai oai \ w c y 



) (84) 



and are oscillating for shorter waves, n > n c (Sakharov' s oscillations). At 
a series of frequencies, the factor sin d is zero, and the growing component 
totally vanishes: no gravitational field perturbations, no time derivatives 
of the perturbations, no energy density perturbations. These modes were 
highly excited, like others, at the end of the i stage, but were stripped off of 
their energy by the very late times of their evolution. In terms of quantum 
mechanics, one can say that these modes have been desqueezed, sent back to 
the vacuum state [23]. [The generating process and associated squeezing 
can be visualized as a depletion of a zero-momentum pump "quantum" 
2uj into a pair of quanta which have the half-frequency uj and oppositely 
directed momenta. Then, the desqueezing can be thought of as the reverse 
process.] The position of zeros is determined by 

-^(m-m) =nk, fc = l,2,3... 

or, approximately, by (n/n c ) = irk. If one defines the distance x traveled by 
sound waves between the barriers at i] = i]i and rj = i]2 by x = a-^(r]2 — rji), 
the zeros arise when x is covered by an integer number of half-waves, 
x = ^k. The first zero in the spectrum of the growing component arises at 
k = 1 which corresponds to the present day scale of the order of 70 Mpc. 

The numerical values of C\ and C 2 as functions of n are controlled by 
the n-dependent function fi(r]i). For simple scale factors (6) and solutions 
(40), n(r]i) is given by Eq. (81) where the value of A\ is determined by 
quantum mechanics, as will be discussed in Sec. VIII. The function n(r]i) 
is exactly the same as the one used for gravitational wave calculations [17]. 
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This allows us to make certain comparisons of density perturbations with 
gravitational waves. For instance, the growing component of gravitational 
waves taken at the beginning of the m stage, h gw (rj2), has the following 
amplitude in the low frequency limit, n <C n c : 



This number is 3\f3ir times larger than C±, Eq. (84), which gives the am- 
plitude h for density perturbations in the same limit. We can also compare 
the growing component of hi with the growing component of gravitational 
waves. Let us take 77 = r]R and n <C nn where nu = ^/(f]R — Vm) corre- 
sponds to the wavelength equal to the Hubble radius at i] = ijr. One can 
derive 



As we see, gravitational waves and density perturbations "enter" the (time- 
dependent) Hubble radius with approximately equal amplitudes, regardless 
of the numerical values of parameters l , (3 describing the i stage. 

According to our definitions, see Sec. VIII, the Fourier component of the 
quantized field includes the factor l p i/y/2n in addition to h(rj) or hi{rj). We 
can give an estimate for the "characteristic" amplitude h(n) ~ nl v \C\ of the 
h field, which is a substitute for a more rigorously defined expectation value 
of the dispersion (square root of the variance) of the field per logarithmic 
frequency interval. Combining Eqs. (84) and (81) we can find in the low 
frequency limit n <C n c : h(n) ~ {l P i/lo)n^ +2 i-e., exactly the same behavior 
as for gravitational waves. The growing components of hi(rj) and Se/e (r]) 
contain the additional factor n 2 (rj — rj m ) 2 which gives two extra powers 
of n in their spectra. There is nothing spectacular about the De Sitter 
case (3 = — 2. The derivative of the Hubble parameter can be arbitrarily 
close to zero at the time when the wavelength of our interest leaves the 
Hubble radius at the i stage. The perturbation will have a finite, not infinite, 
amplitude today. 

7. Density and Rotational Perturbations 
in the High-Frequency Limit 

The normalization of the perturbations is determined by quantum me- 
chanics. We intend to amplify the zero-point quantum fluctuations of the 
primeval matter which is, in our case, the scalar field (1). Before the am- 
plification, the frequencies of the fluctuations were much higher than the 
frequency of the gravitational pump field. To the modes of our interest the 
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surrounding space-time seemed at the beginning almost flat. As a prepara- 
tion for quantization, we will first consider density and rotational pertur- 
bations in matter placed in the Minkowski space-time, that is when gravity 
is totally neglected (a(rj) = 1 in Eq. (2)). 

The deformation of an elastic medium is usually described [24] with the 
help of a displacement three- vector Uj(t,x) which can be written as 

«i = Z(t)Q,i + &(t)Qi . (85) 

The scalar function Q is defined by Eq. (8), the vector function Qi is defined 
by the equations 

Qi,k' k + n 2 Qi=0, Q\i = . (86) 
The deformation tensor is 

u ik = ^{u itk + u k}i ) = ^(Q,i,k + Q,k,i) + ^0(Qi,k + Qk,i) 

and its trace is u = u % ^ = —^n 2 Q. The stress tensor can be written in the 
general form 

Oik = s{Q,i, k + Q,k,i) + PCtO(Qi,k + Qk,i) + n 2 {2s - pc 2 £)Q5 ik (87) 

where p is density, q is the longitudinal velocity of sound, q is the trans- 
verse (torsional) velocity of sound, and s is arbitrary function of time. The 
equations of motion 

d 2 Ui _ do\ 
^ dt 2 dx k 

reduce to the oscillatory equations for elastic waves. In terms of 77-time, 
they can be written as 

e"+n 2 j£ = 0, 0» + n *$e = . (88) 

We shall now relate the theory of elasticity with the theory of cos- 
mological perturbations. We shall consider the high-frequency limit of the 
perturbations, that is when the scale factor a(rj) is almost constant and its 
variability can be neglected in comparison with frequencies of the waves. As 
we know, the perturbed components of the energy-momentum tensor for 
density and rotational perturbations have the general form (see Section II 
and Ref. [2]): 
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(89) 

The differential conservation laws T@ * = can be reduced in the high- 
frequency limit to the following equations: the a = component gives 
the equation — e[ + £'n 2 = 0, which integrates to 

n 2 e = ei , (90) 

the a = i components give 

e"+ Pl = 0, 9" + X = ■ (91) 

Equations (90) and (91) can, of course, be obtained from the perturbed 
Einstein equations as well. 

The stress tensor is connected with the perturbed components T k 
by (Ji k = —pc 2 Ti k (a(rj) = 1). From the comparison of Eqs. (87) and (89) 
one finds 

Pl = n 2 |e, X = n 2 p, s = -pc 2 ^. (92) 

With these expressions for p\, x, Eqs. (91) coincide with the wave equations 
(88) . In cosmology and theory of elasticity, we are dealing essentially with 
the same physics. 

The quantization of density and rotational perturbations should be 
based on Eqs. (88). Rotational perturbations have been considered else- 
where [2]. The quantum-mechanically generated rotational perturbations 
can contribute to the CMBR anisotropy and may be important for the 
smaller scale astrophysics. (One should be aware, though, that there exists 
also an alternative view on the subject according to which rotational per- 
turbations are "irrelevant for cosmology" [25].) We will concentrate here 
on density perturbations. For simple models of matter, such as perfect fluids 
and scalar fields, the functions x an d Pi vanish. One is left with the single 
variable £ and isotropic stresses. The quantization of these oscillations in 
an elastic material placed in the Minkowski space-time would lead to the 
notion of phonons. 

There is no wonder that in the case of scalar field matter the role of £ is 
played by ip±. If one writes £(r/) = £ e~ mr > and ^1(77) = 9910 e~ mv , Eq. (24) 
gives 

m£ = VoPio ■ (93) 

In the high-frequency limit (large n), the first term in Eqs. (23) and (25) 
dominates. With the help of Eq. (93) one derives e\ = pi = n 2 £. In other 



1 Pi 

a 2 2n 2 



Q/ + Q' 



+ ^(pi+Pi)Si k Q. 
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words, for the high-frequency scalar field perturbations, the velocity of 
sound is almost equal to the velocity of light. The perturbations behave 
as massless scalar particles. 

It is the scalar field oscillations that should be normalized by ascribing 
a "half of the quantum" to each mode. Due to the Einstein equations the 
scalar field perturbations are accompanied by the gravitational field per- 
turbations. In the high-frequency limit, the second term in Eq. (43) can 
be neglected, the normalization of tp\ transfers to the gravitational field 
variable fx and ultimately to h and hi. This is how we will know the initial 
amplitude for density perturbations. 



8. Quantization of Density Perturbations 



In the limit of a free massless scalar field placed in the Minkowski space- 
time, we have for each mode: 

<Vk = fiQ + tiQ* = Mt) e lky + <pl(t) e~ iky . 
The total field can be written as 



6<p(t,y)=C- 



1 



OO 1 

d 3 k^L= 

-oo V^Wk 



c k e~ iWkt e iky + 4 e iWkt e~ iky l . (94) 



(2tt) 3 / 2 

The normalization constant C is to be found from the requirement 



1 



<0| / eo! 3 y|0) = Jft / ^ku; fe (0| Ck 4 + 4c k |0) , 

J —oc ^ J —oc 

where e is the energy density of the field. Since in our case 

L Wo)' + (<V,l) 2 + (<V,2) 2 + (<V, 3 ) 2 

we derive C = cVh. Obviously, the normalization coefficient C includes the 
Planck constant h but does not include the gravitational constant. 

To write the field operator (94) in the curved space-time (2) one should 
make the following replacements: 



y = a(r/)x, k 



a(r?) 



n, w k 



cn 
a(r?) 



. c k = Hv)? /2 Cn 



The field operator takes on the form 



<M»7> X ) 



1 



afa) (2vr) 3 / 2 



OO 1 

d 3 n^[c n e~ inv e inx + t£ e inr? e~ inx ] . 
-oo v2n 



(95) 
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This expression is only valid in the high-frequency limit, when the field 
can be regarded as free and its non-adiabatic interaction with gravity can 
be neglected. When the interaction becomes important, the time depen- 
dence of the field ceases to be so simple. The operator c n e~ mv , should be 
replaced by c n (77) , and the evolution of c n (77) , Cn (77) , should be found from 
the Heisenberg equations of motion. 

The matter perturbations are accompanied by the gravitational field 
perturbations. Due to the Einstein equations they are linked together and 
form, in a sense, a united entity. As we have seen in Sec. Ill, the entire 
dynamical problem at the i stage reduces to a single wave equation (38) for a 
single variable n{rj). All perturbations can be found from a given solution to 
this equation. It follows from Eq. (43) that ix(rj) « \/2k aipi(r)) in the high- 
frequency limit. Since_the positive frequency scalar field n-mode solution 
is ipi(r)) ~ (l/a(r]))Vchc n e~ tnv this leads to /i(r/) ~ A^/ti l p ic n e~ mv . Note 
that the normalization coefficient includes the gravitational constant and 
is proportional to the Planck length l p i = (Gfr/c 3 ) 1 / 2 . 

Having in mind the basic wave equation (38), we can now introduce 
the "fundamental" scalar field $(77, x) which describes the whole quantum 
system interacting with the gravitational pump field: 



1 f°° 



c n (r7)e inx + C t n (7 ? )e- inx 



• (96) 



The annihilation and creation operators c n (r]), Cn(rj) are governed by the 
Heisenberg equations of motion 

—j— = -i[cn(n), H] , — - — = -iK{7]), H] . (97) 

The interaction Hamiltonian H is given by 

H = nc^Cn + racL n c_ n + 2a(7])clci n + 2a* (r/)c n c_ n , (98) 
where the coupling function a{j]) is 



(7(77) 



2 0^7 



Equation (98) demonstrates explicitly the underlying parametric interac- 
tion of the field oscillators with the pump field. This is simply a generaliza- 
tion of a theory previously developed for gravitational waves (for a review, 
see Ref. 19). 
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The common way of solving Eqs. (97) and (98) is to write the operators 
in the form 

c„(r?) = u n (r))c n (0) + v n (r])c j _ n (0) , 4 fa) = < (77)4(0) + <(r?)c_ n (0) , 

(99) 

where c n (0), 4(0) are the initial values of the operators taken at some 
77 = 770 long before the interaction became effective, and [c n (0), 4i(0)] = 
<5 3 (n — m) The classical complex functions u n (r)), v n (rj) (do not mix up 
with the functions u(rj), v(rj) introduced in Sec. Ill) obey the condition 
\un\ 2 — \v n \ 2 = 1 and satisfy the equations 

iu' n = nu n + 2av* n , iv' n = nv n + 2<ju* n (100) 

with the initial data u n (0) = 1, 7j n (0) = 0. If one introduces fi n (rj) = 
Un(v) + v n(v)i it follows from Eqs. (100) that the function n n {rj) should 
satisfy precisely the Eq. (38). The initial conditions for n n (rj) in the high- 
frequency limit |n?7| — > 00 are fJ-niv) ~^ e _i "^ ??_??0 \ n' n (f]) — > —ine~ m ( ri ~ ri0 \ 

For each mode n there exists the vacuum state |0 n ) defined by the 
condition c n (0)|0 n ) = 0. As a result of the Schrodinger evolution with 
the Hamiltonian (98), the initial vacuum state |0 nj _ n ) = |0 n )|0_ n ) trans- 
forms into a multiparticle two-mode squeezed vacuum state (see Ref. 19 
and references cited therein). In other words, the perturbations (waves) are 
generated in correlated pairs. Squeezing is a phenomenon, not a formalism. 

At this point, it may be useful to discuss the difference between squeezed 
states and other quantum states. It is sometimes said that the strongly 
squeezed states are "almost classical" since their wave functions can be 
written in the WKB form and the noncommutativity of the creation and 
annihilation operators is not important. The implication seems to be that 
there is no need for quantum-mechanical treatment of cosmological pertur- 
bations, everything can be done "classically". I think that these arguments 
are misleading. It is true that the strongly squeezed quantum states are the 
multiparticle states, the mean number of quanta in these states is much 
larger than one. Long ago, but not now, this was considered as a criterium 
for the state to be "classical". However, there are many quantum states 
which are the "WKB classical" and whose occupation numbers are very 
big and, yet, the statistical properties of these states can be entirely dif- 
ferent. Two quantum states may have exactly the same, and big, mean 
number of quanta but, along with this, they may have drastically differ- 
ent variances of this number. This is precisely what is true with regard to 
squeezed states, as compared with the "most classical" quantum states - 
coherent states. 
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Exercise j. 

Show that the mean value of the number operator N = a) a in a squeezed 
vacuum state is (N) = sh?r. Let the squeeze parameter r be r > 1. Con- 
sider a coherent state with exactly the same mean number (N). Show that 
in the coherent state the variance {(AN) 2 ) is ((AiV) 2 ) = (N) while in the 
squeezed vacuum state it is ((AN) 2 ) = ^sh 2 2r « (N) 2 . 



The statistical properties of the field are determined by c n (rj), c a (r]). 
By using Eq. (99) one can rewrite Eq. (96) in the form 

c n (0K(?/K nx + ct(oK(77) e - in x 

(101) 

where the functions /x n (^/), lAiiv) should be taken with the appropriate ini- 
tial conditions discussed above. For simple solutions (40), the initial condi- 
tions translate into the requirements 

A x = * /| e <(™»+¥), A 2 = iA 1 e- i *P . (102) 

cos fin V 2 

The quantized gravitational field perturbations are expressible entirely 
in terms of the $(r/, x) field (96). There are many components of hij but 
there is only one sort of creation and annihilation operators. Let us intro- 
duce new notations 

l 2 

h(v) = h(r]), h(v)=hi(r}) 

and the polarization tensors 

1 2 „ 1 2»J 

Pij = $ij , Pij = -nirij/n , and PijP = -1 . 
The field operator hij(r],-x.) can be written as 

c n (0)Le inx + 4(0)£e- inx 

(103) 

where the classical complex functions hn(v) should be derived from /J, n (rj) 
through the equations and initial conditions already discussed. A similar 
expression can be written for the operator of the energy density perturba- 
tion Se/e . Equation (103) is the starting point for the calculation of the 
expected angular anisotropy in CMBR. 



$(77, x) = Ay/nl 



P l 



(2vr) 3 / 2 



d 6 n 



\ f 00 1 s 

4 ^"(2rf7I L^TTE ?>«<»> 
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9. Variations of the CMBR Temperature Caused 

by Density Perturbations of Quantum-Mechanical Origin 

The photons of CMBR are emitted at r] = tje and are received by us at 
r] = r]R. A particular direction of observations is characterized by the unit 
vector e k = (sin 9 cos (/), sin 9 sin (ft, cos 9). In absence of perturbations, tem- 
perature of CMBR seen in all directions would be the same, T. Gravi- 
tational field hij associated with the density perturbations at the m stage 
causes a variation of the temperature with respect to the unperturbed value 
T [26]: 

dhi 



(104) 



where w\ = t]r — t]e, and dhij/drj is taken in the comoving synchronous 
coordinate system along the path x k = e k w, i] = r]R — w. In case of small 
perturbations, which we are actually dealing with, the emission time t)e can 
be regarded as being one and the same for all directions. Since the scale 
factor satisfies the approximate relationship a{r)E) / a{r)R) ~ 10~ 3 , and the 
rj time can be chosen in such a way that t]ji — r] m = 1, the quantity w\ is 
close to 1. We will use w± = 1. The wavelength equal to the present day 
Hubble radius In corresponds to n# = Ait. 

For the quantized perturbations, the temperature variation 5T/T, 
Eq. (104), becomes a quantum- mechanical operator. Since it is only the 
longitudinal-longitudinal part of hij that participates in producing 5T/T, 
we can write 



5T 



( e *) = -2^l pl 



j^wl dw L 



d 3 n 



me 



i\2 



Cn(fyfn(VR ~ W)e 



:!!kekw + ci(0)f:( VR -w)e- m ^ kw 



(105) 



where 



fn(VR 



w 



1 dh^ij) 



In di] 



The individual observed distributions of the CMBR temperature over 
the sky should be compared with theoretical predictions based on the 
quantum-mechanical expectation values. The mean value of 5T/T(e k ) is 
obviously zero, (0\5T/T (e k )\0) = 0. The variance (0\8T/T (e k )5T/T(e k )\0) 
is not zero but does not depend on the point and direction of observations. 
To study the angular distribution of the temperature variations one should 
construct the angular correlation function K for two different directions e k 
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and e 2 , efe 2 5ki = cos 5: 

KHQ|f^)f^)|o) . 

By manipulating with the product of two expressions (105), one can derive 

1 f'^ /"OO /'TT 

if = 4-7r/p Z - — / / cfe / n 2 \f n \ 2 dn / sin^c?^ 
(27r)' 3 Jo Jo Jo Jo 

* ( " ,e ' l)2 ( ™ ,e ' !)2 co S (, H ,C*)# , (106) 



/ 

Jo 



n 2 n 2 



where C k = e i w ~ e 2 w. A lengthy calculation of the integrals over the 
angular variables </>, 6 gives the following result: 



[\m9de f 2W ^ e p (gjfp cos(n fc C fc )# 
Jo jo n 2 

4vry||cos 2 ( 5(nC)- 1 / 2 J 1/2 K) 

+ (l-5cos 2 ,5)(nC)- 3/2 J3/ 2 K) 

+ 2 cos <5(1 - cos 2 <5) (ot) (nw) {n()~ 5/2 J 5/2 (n() 

+ 4(3cos 2 ,5-l)(nC)- 5/2 J5/2K) 

+ 8 cos 5(cos 2 5- l)(nw)(nw)(n()~ 7/2 J 7/2 (n(;) 

+ (cos 2 5-l) 2 (nw) 2 (nw) 2 (n()- 9/2 J 9/2 (n() } (107) 



where C = ("W 2 — 2unt> cos 5 + w 2 ) 1 ^ 2 . For further calculations one may use 
the following formula (valid for half-integer v): 

( j i a\ /^TX^t , ; ,Ju+k(nw) J u+k {nw) d v ~ x l 2 

(n() u Jv{nQ = V2vr ^ (y + k)— , P k+u _ 1/2 (z) 

(nw) u (nw) u dz u 1 ' z ' 

(108) 

where z = cos 5 and Pi(z) are the Legendre polynomials. (I derived and 
used this formula in course of studying the gravitational wave [17] and ro- 
tational [2] perturbations, but I believe that this formula may exist some- 
where in the previously published literature.) With the help of Eq. (108), 
one can rearrange the correlation function K to the following final expres- 
sion: 

oo 

K = l 2 pl J2KiPi(cos5) , (109) 
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where 



Ki = (21 + 1) / \I ln \'- 
Jo 



dn 



I, 



In 



n fn{m ~ x/n) 

,1/2 



1(1-1) 



1 



Ji +l (x) + -Js +l (x)Ux, (110) 



and x = nw. The derived formula for K is general and can be used with 
arbitrary function /„. In practice, the limits of integration over n are deter- 
mined by the frequency interval within which the perturbations were really 
generated. 

As we see, the decomposition of K consists of all multipoles including 
the monopole, I = 0, and dipole, / = 1, terms. To carry out the calculations 
up to a concrete number, we will consider scale factors (6) and solutions 
(40). 

As we know, the growing and decaying components of hi(rj) are of the 
same order of magnitude at n = r\i- However, the decaying component is 
decreasing since then and can be neglected in the calculation of 5T/T. For 
the coefficient C\ responsible for the growing solution, see Eq. (82), we have 



|d| 2 



^(/?)|V^sin 2 ™ 



(111) 



where 

\m\ 2 



2^+5 cos PttT(P + ^) 



|V(/3)| 2 = 1 for/? = -2 



The expression for K\ takes on the form 
1 



Ki = 



100Z 2 , 1 



W)\ 2 n 2 c (2l + l)j n 



2(3+1 







n 
n c 



sin —(tin) dn. 



(112) 



where 



r n n — x 
= Jo 1^ 



1(1-1) 



- 1 



Ji +l (x) H — Jx +[ (x) >dx . (H3) 

2 X 2 | 



We will now estimate the contribution of long waves, n < 1, to the lower 
order multipoles K%. The integrals i\ n should be calculated separately for 
I = 0, / = 1, and / > 2, because of the factor /(/ - 1) in Eq. (113). For / > 2, 
the term with this factor dominates. The approximate expression for (i; n ) 2 , 
I > 2, takes on the form 



(iln) ~ A 



(114) 
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where 




For / = and I = 1, the term with the factor 1(1 — 1) does not contribute 
to (iin) 2 ''■ The result still have the form of Eq. (114) but n 21 should be 
replaced by j^n 6 for I = 1, and by ^n 4 for I = 0. These results are in 
full agreement with Ref. 13: in the limit of long waves, the monopole and 
dipole contributions of an individual wave are suppressed; the monopole 
component is of the same order of magnitude as the quadrupole component, 
while the dipole component is further suppressed by an extra power of n. 

We should now use (ii n ) 2 for the calculation of K\. This is where the 
spectrum of the perturbations comes into play. We can write for I > 2: 

K l = T ^\m\ 2 m + ^ I'n^^dn , (115) 

and the appropriate replacements discussed above should be made for / = 1, 
/ = 0. Since we are working in the limit of small n, the integration over n 
cannot be extended to the values n > 1. However, typically, short waves 
contribute little to the lower index multipoles. The values of n up to n « 21 
are more important, but for the purposes of simple evaluation we restrict 
the integration by n = 1. 

The suppression of the monopole contributions of individual waves saves 
us from a big trouble. If Eq. (115) were true for I = 0, the monopole term 
Kq would be power-law divergent in the limit of n — ► for all (3 < — 2. 
In order not to be in conflict with the finite observed 2.7 K temperature, 
we would need to resort to the fine tuned minimally sufficient duration of 
the i stage, in which case the long waves with this spectrum are simply 
not being generated. From the correct Eq. (115) it follows that the danger 
of divergence for Kq and K2 arises only in models with (3 < —4. (The 
quadrupole anisotropy produced by gravitational waves does also diverge 
in these models [17].) Thus, the interval — 2 > [3 > — 4 is potentially 
allowed. 

By the way, the level of the quadrupole anisotropy puts sensitive quali- 
tative constraints on the conjectures [31] about the "very inhomogeneous" 
Universe on scales "far beyond the visible horizon" . Even small linearized 
perturbations can not have totally unrestricted amplitudes and spectra on 
scales longer than Ih, see [13,32,33]. 

We will now introduce the notations I p i\/Kq = M, l v i\[K\ = D, 
lpiVT^2 = Qi and will compare M, D, and Q. For the quadrupole Q, one 
can find from Eq. (115) 

Q ^'t Hm 7m ■ (116) 
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The monopole M and dipole D are related with Q by 

M _ \/5 L> _ /T / yg + 4 
Q"~~2~' Q ~ V20V/3 + 5 ' 

We do not have an observational access to the unperturbed temperature 
T, but whatever is the measured Q, we can expect that a correction of 
about the same magnitude as Q is included in the measured T. The same 
is true for the dipole component D (there is little doubt, however, that the 
overwhelming part of the measured dipole anisotropy is accounted for by 
our peculiar motion). 

We will now compare the contributions of density perturbations and 
gravitational waves to the components K\ of the correlation function K 
in the long wave limit, n < 1. We should compare Eq. (115) with the 
analogous expression for gravitational waves [17]. The ratio of the gravity 

g d 
wave contribution Ki to the density contribution Ki has the form 

ki (i + m + 2)(2i + if 




21(1 - 1) 



We see that the ratio is independent of the parameters l Q , (3 describing the 

9 d 

i stage. For the quadrupoles Q and Q, we have in the long wave limit 
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that is a somewhat larger contribution of gravitational waves. It is necessary 
to take into account also the shorter waves in order to get a more accurate 
estimate. 

In conclusion, there is no dimensionless ratios that could be adjusted in 
such a way that the contribution of density perturbations to the quadrupole 
anisotropy would be much larger than the contribution of gravitational 
waves. These contributions are of the same order of magnitude while nu- 
merical coefficients are somewhat in favour of gravitational waves. At the 
same time, the very generation of density perturbations (and rotational per- 
turbations) is more problematic than the generation of gravitational waves. 
On these grounds one can conclude that if the observed large-angular-scale 
anisotropy of CMBR is caused by cosmological perturbations of quantum- 
mechanical origin (what else?), they are, most likely, gravitational waves. 
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